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ABSTRACT 

We construct merger trees based on the extended Press-Schechter theory (EPS) in 
order to study the merger rates of dark matter haloes over a range of present day 
mass (1O 1O M < M < 1O 15 M ), progenitor mass (5 x 1CT 3 < £ < 1) and redshift 
(0 ^ z ^ 3). We used the first crossing distribution of a moving barrier of the form 
B(S,z) = p(z) + q(z)S 7 , proposed by Sheth & Tormen, to take into account the 
ellipsoidal nature of collapse. We find that the mean merger rate per halo B m /n 
depends on the halo mass M as M - 2 and on the redshift as (dS c (z)/dz) 11 . Our results 
are in agreement with the predictions of N-body simulations and this shows the ability 
of merger-trees based on EPS theory to follow with a satisfactory agreement the results 
of N-body simulations and the evolution of structures in a hierarchical Universe. 
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1 INTRODUCTION 

Although numerical experiments are the most powerful 
methods to study the formation of structures, the devel- 
opment of analytical or semi-numerical methods is very im- 
portant, since they help to improve our understanding of the 
physical processes during the formation. 
A class of analytical methods is that based on the ideas of 
Press & Schechter (1974) and on their extensions (Bond et 
al. 1991, Lacey & Cole 1993): 

The linear overdensity computed at a given point of an ini- 
tial snapshot of the Universe fluctuates when the smooth- 
ing scale decreases. This fluctuation is a Markovian pro- 
cess when the smoothing is performed using a top-hat win- 
dow in Fourier space. For any value of the smoothing scale 
R the overdensity field is assumed to be Gaussian. The 
mass M contained in a given scale R depends on the win- 
dow function used. For a top-hat window the relation is: 



M 



4 „ 



R'- 



-R 3 , where pb,i and fi m ,i are the 



values of the mean density and the density parameter of 
the Universe, G is the gravitational constant and Hi is the 
Hubble's constant. The index i indicates that all the above 
values are calculated at the initial snapshot. Mass disper- 
sion a 2 at scale R is a function of mass M and is usually 
denoted by S, that is S(M) = a 2 [R(M)]. Let the random 
walk of the overdensity cross for first time a given barrier 
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B(S,z) at some value So of 5. The mass element associated 
with the random walk is considered to belong to a halo of 
mass Mo = S~ 1 (So) at the epoch with redshift z. The dis- 
tribution of haloes, at some epoch z, is related to the first 
crossing distributions, by the random walks, of the barrier 
that corresponds to epoch z. 

The simplest form of the barrier comes from the spherical 
model. It is well known that in an Einstein-de Sitter Uni- 
verse, a spherical overdensity collapses at z if the linear ex- 
trapolation of its value up to the present exceeds S 3C « 1.686, 
and this value provides a first reasonable choice for the bar- 
rier. The involved quantities (density overdensities, disper- 
sion) are usually extrapolated to the present epoch and thus 
the barrier in the spherical collapse model is written in the 
form B(S, z) — 1.686/D(z), where D(z) is the growth factor 
derived by the linear theory with D(z = 0) = 1. The form 
of the spherical barrier permits the analytical evaluation of 
the first crossing distribution f(S). 

Despite the simplicity of the spherical model, its results 
agree relatively well with the results of N-body simulations 
(e.g. White, Efstathiou & Frenk 1993; Lacey & Cole 1994; 
Gelb & Bertschinger 1994; Bond & Myers 1996). Deviations 
appear in the resulting mas s func tions at both high and low 
masses. Sheth & Tormen (1999) consider a barrier of the 
form B(S,z) = p(z) + q(z)S"< with p{z) = 0.8405 c (z), q(z) = 



0.505<5 c (z) 



where 8 c {z) = 1.686/D(z) and 7 = 0.615, 



in order to describe the ellipsoidal form of collapse. The first 
crossing distribution f(S, z) that results from the ellipsoidal 
model works better than the spherical one. For example, 
Yahagi et al. (2004) showed that the multiplicity function 
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resulting from N-body simulations is far from the predic- 
tions of spherical model while it shows an excellent agree- 
ment with the re sults of the ellipsoidal model. On the other 
hand, Lin et al. (2003) compared the distribution of forma- 
tion times of haloes formed in N-body simulations with the 
formation times of haloes formed in terms of the spherical 
collapse model. They found that N-body simu lation s give 
smaller formation times. Hiotelis & del Popolo (|2006h used 
merger trees to show that the ellipsoidal collapse model leads 
to formation times that are shifted to smaller values relative 
to the spherical collapse model. Thus, a better agreement 
with the predictions of N-body simulations is achieved. The 
distributions of formation times are studied in more detail 
by other authors (Giocoli et al. 2007). We note that in the 
cases of a barrier with 7 = 1, or 7 = 0. 5, the correspond- 
ing first crossing distribution can be found analytically (e.g. 
Mahmood & Raje sh 200 5). For the above non-linear barrier, 
Sheth & Tormen (J2002J) proposed an analytical expression. 
It is shown that this analytical expression approximates well 
the exact expression found as a numerical solution of an in- 
tegral equation (Zhang & Hui 2008 ). 

The constrained first crossing distribution is given by the 
relation: 



f(S,z/S, z )dS- 



1 \T(S,z/S ,z )\ 
2lF (AS) 3 / 2 
(AS) 2 



cxp 



2AS 



d5 



where AS = S — So, AB = B(S,z) 
function T is given by: 

T(S, z/So, zo) = B(S, z) — B(So, zo) + 



(1) 



B(So, zo) and the 
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Given that a mass element is a part of a halo of mass Mo 
at redshift zo the probability that at higher redshift z this 
mass element was a part of a smaller halo M is given by 
Eq.(l). The unconstrained expression, f(S, z), results by set- 
ting So = B(So,zo) = 0. The analytical expression is very 
useful since it allows the construction of merger trees. 
The purpose of this paper is: 

(i) To construct merger trees able to give merger rates of 
dark matter haloes. 

(ii) To compare these merger rates with those predicted 
by the results of N-body simulations. 

(iii) To extend the calculations to scales that are not ac- 
cessible by numerical simulations and finally, 

(iv) to study the role of some of the main parameters 
involved. 

In Sect. 2 we give a brief description of the tree code used. 
Then, the definition of merger rates is presented and the 
analytical formulae, predicted by other authors from the re- 
sults of N-body simulations, are given. 
In Sect. 3 we present our results. 



2 TREE CONSTRUCTION, DEFINITIONS OF 
HALO MERGER RATES FITTING 
FORMULAE 

2.1 Tree construction 

Merger-trees used in this paper are constructed using 
Eq.(l). Let us assume a descendant halo of mass Md at 
redshift Zd- We study its past by the following procedure: 
A new larger redshift z v is chosen. This is done by solving 
for z p the equation S c (z p ) — S c (zd) = D, where D is a 
constant (one of the parameters of the algorithm). Then, 
a value S p is chosen from the distribution (1). The mass 
of the progenitor is M p = S~ 1 (S P ). This progenitor is 
accepted if its mass is larger than a lower limit M m ;„ and 
smaller than the mass left to be resolved. If M p is less than 
M min then M p is added to a sum that is named M accr . 
The mass left to be resolved is, at the choice of the k-th 
progenitor, M; e /t = M d — Xa=i _1 M p, 1 ~ M ».ccr- If M ie f t 
is larger than M m i„ we proceed to the selection of a next 
progenitor, while if M; e / t is smaller than M m in we proceed 
to the next redshift. It is obvious that in such construction, 
the number of progenitors can be larger t han tw o, despite 
the original assumption of Lacey & Cole (1993), and this 
leads to a better representation of the distribution of 
progenitors. Our algorithm is ba sed on the "N-branch" 
idea of Somerville & Kollat (|l999h . but extended in order 
to incorporate aspects of the ellipsoidal collapse results. A 
complete description of the t ree construction is given in 
Hiotelis & del Popolo (2006). The comparisons between 
the predictions of the tree and analytical predictions of the 
distribution of the number of progenitors show that the tree 
method is reliable in following the evolution of structures. 
Various tree-construction algorithms have been presented 
in the literature (Cole 1991; Kauffmann & White 1993; 
Sheth & Lemson 1999; Cole et al. 2000; Neinstein & Dekel 
2008). The accuracy of these algorithms is usually against 
simpli city. For example, the algorithm of Neinstein & Dekel 
(2008) requires the solution of several differential equations 
with nontrivial boundary conditions. 



2.2 Definition of merger-rates 

We examine one descendant halo from a sample of Nd haloes 
with masses in the range Md, Md + dMd present at red- 
shift Zd- For a single halo the procedure is as follows: Let 
M p ,i, M Pt 2---M Pik be the masses of its fc progenitors at red- 
shift z p > Zd- For matter of simplicity we assume that the 
most massive progenitor is M Pl i. We define £i = M Pt i/M P: i 
for i ^ 2 and we assume that the descendant halo is formed 
by the following procedure: During the interval dz — z p — Zd 
every one of the progenitors with i ^ 2 merge with the most 
massive progenitor i = 1 and form the descendant halo we 
examine. We repeat the above procedure for all haloes in the 
range Md, Md + dMd found in a volume V of the Universe. 
Then, we find the number denoted by N of all progenitors 
with £i, i ^ 2 in the range (£, £ + d£) and we calculate the 
ratio N/(VdzdMdd£). We define the merger rate B m as fol- 
lows: 



B m {M d ,i,z p : z d ) 



N 



VdzdM d d£ 



(3) 
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Let the number density of haloes with masses in the range 
M d ,M d + dM d at z d be n(M d ,z d ) = Nd j^ d) ■ The ratio 
B m /n = N/(N d dzd£) measures the mean number of merg- 
ers per halo, per unit redshift, for descendant haloes in the 
range M d , M d + dM d with progenitor mass ratio £. 
We note that the definition of mean merger rate is exactly 
the same as in Fakhouri & Ma (|2008l ) (FM08 hereafter). We 
also use the assumption that all progenitors merge with the 
most massive one (see FM08 for a discussion of this assump- 
tion.) 

Lacey & Cole (|l993h showed that in the spherical model the 
transition rate is given by: 



r (M — ♦ M d /z d )dM d = (|) ] 
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dM d 



(4) 



This provides the fraction of the mass belonging to haloes 
of mass M that merge instantaneously to form haloes 
of mass in the range M d ,M d + dM d at z d . The prod- 
uct r ■ f sc (M,z d )dM, where f sc (M,z) is the unconditional 
first crossing distribution for the spherical model, gives the 
above fraction of mass as a fraction of the total mass of 
the Universe and successively multiplying by (pb/M) ■ V 
the number of those haloes is found. Then, by dividing by 
(pb/M d ) ■ V ■ f S c(M d , z d )dM d (that equals to the number of 
the descendant haloes) we find: 



N 
N d dz 



2 M d 



da(M) 
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(5) 



Assuming a binary merge, where £ is the ratio of the small 
progenitor to the large one (£ = (M d — M)/M), using dM = 
and substituting in (5) we have the final expression 
for the binary spherical case, that is: 



Br, 



N 



N d dzd£ 



M da(M) \dS c (z) 



TV a 2 (M) dM 



dz 



a 2 (M d 



— z d 
-i -3/2 



<T 2 (M) 



(6) 



2.3 Fitting formulae 

FM08 analyzed t he res ults of the Mille nnium simulation 
of Springel et al. (I2005D . Stewart et al. (|2008l ). (SBBW08 
hereafter), used their high- resolution N-body simulations 
to study the merger rates of dark matter haloes. Fitting 
formulae proposed by the above authors are separable in 
the three variables, mass M d , progenitor ratio £ and redshift 
z. These formulae are of the form: 



B { M d ,$,z p : z d ) 
n(M d ,z) 



A- F{M d )G(i)H{z) 



(7) 



FM08 proposed A 

i a 31 



C 2 exp 



(!)' 



0.0289, F(M d ) = 



,H(z) 



where the values 



of 



z = z d 

-,12 , 



dz 

the parameters are M = 1.2 X W L2 M Q ,A = 0.0289, £ 
0.098, ai = 0.083, a 2 = -2.01, a 3 = 0.409, a 4 = 0.371. 
On the other hand, SBBW08 proposed: A = 0.27, F{M d ) = 

W^G(0 = (l-C) b3 - 1 [(b3-&2)C + b2]/C i,2+1 ,^W = 
(¥i) b4 _ where M = 1O 12 ^ _1 M ,6i = 0.15,62 = 0.5, b 3 = 
1.3, b 4 ~=2. 

Note that the formulae proposed by the above authors show 
some significant differences. First, the dependence on red- 
shift z differs in the above two formulae. Although the quan- 
tity dS c /dz does not vary significantly with redshift, the ex- 
ponents CI4 and 64 can cause significant differences in the 
merger rates. Second, exponents ai and 61 that define the 
dependence on the mass are quite different. In the approxi- 
mation of FM08 mean merger rates are practically indepen- 
dent on mass. 

Additionally, the above formulae show significant differences 
in their slope at small and large values of £. The logarithmic 
slope of G, dlnG(£)/dln£, at £ = is a 2 = -2.01 for the 
FM08 model, and -b 2 - 1 = -1.5 for the SBBW08 model. 
For £ — > 1 the above logarithmic slope is a 2 + ci3 ( 1 /f )" 3 — 
—0.952 for the FM08 model and it tends to —00 in the for- 
mula of SBBW08. 

Therefore, it is interesting to study merger rates predicted 
by merger-trees and to compare their characteristics with 
those of the above fitting formulae. 



3 RESULTS 

We used a flat model for the Universe with Q m ,o = 0.3 and 
Qa.o = 0.7 and a power spectrum proposed by Smith et al. 
(|l998T ) given by: 



P(k) 



Ak n 



[1 + aik 1 / 2 + a 2 k + a 3 k 3 / 2 + ai k 2 ] b 



(8) 



The values for the parameters are: n = 1, ai = 
-1.5598, a 2 = 47.986, a 3 = 117.77, a 4 = 321.92 and 
b — 1.8606. Smoothed fields are calculated using the top- 
hat window function. The constant A of proportionality, is 
found using the procedure of normalization. We used two 
different values for the normalization namely as = o(R = 
8/i _1 Mpc) = 0.9 and 1 respectively. We also use a sys- 
tem of units where M unit = 10 12 h' 1 M Q , R unit = /i _1 Mpc 
and iunit = 1.515 x 10 7 ft -1 years. In this system of units, 
Ho/Humt — 1.5276. We performed a large number of tree 
realizations in order to study merger rates. First, we found 
that good fits are achieved by both of the above formulae but 
for different values of the parameters than those proposed 
by the above authors. We found that merger rates depend 
on the mass of the descendant halo as M d 2 and on the red- 
shift as [dSc/dz] 1 ' 1 . So, in the comparisons that follow, the 
formulae of Eq. 7 are used with ai = b\ = 0.2, a 4 = 64 = 1.1. 
Additionally, in the formula of SBBW08, we use the value 
of 62 = 0.7 instead of b 2 = 0.5. 

In Fig.l, comparisons between the predictions of merger 
trees and those of formulae given by Eq. 7 are shown. The 
power spectrum used is that given by Eq. 8 for as — 1. De- 
tails are given in the caption of the figure. We note a very 
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Figure 1. From top to bottom, dots correspond to the mean 
present merger rates (z = 0) for haloes with present day masses 
lQ 1B h- l M Q , IO^H-^Mq, 10 12 /i- 1 M q and 10 10 h- 1 M Q respec- 
tively, predicted by the formula proposed by FM08 for a\ =0.2 
and 04 = 1.1. Dashed lines are predicted by the formula proposed 
by SBBW08 for 61 = 0.2 and b 2 = 0.7. Solid lines arc the predic- 
tions of merger trees described in the text by Eq.(3) for z d = 
and z p = 0.0556. We used a sample of 10000 present-day haloes 
for each of the three cases and we evolved the system by a single 
time-step for D = 0.05. The minimum mass M m in in every case 
equals to 0.005 times M d . 



good agreement but also a rapid fall of the predictions of 
merger trees for small values of £. We will show below that 
this is a matter of resolution that depends on the value of 

Mmin . 

The predictions of this figure are just after one time step 
of the tree algorithm. For a higher redshift z\ one has to 
use one of the two alternatives: To start with a sample of 
haloes at zi and after a single time step to move to a new 
redshift z p i, to calculate the merger rates, or to start with a 
sample of haloes at the present epoch z = and to calculate 
merger rates after a number of time steps when the redshift 
has a desirable value. Obviously, the second approach tests 
the ability of the merger tree to follow the time evolution 
of structures and it is close to the nature of N-body simula- 
tions. This approach is followed in our calculations. 
Fig. 2 shows merger rates predicted by a merger tree, started 
at z — 0, after 22 time steps, for D = 0.05 and cr$ = 1. Red- 
shifts are z d — 0.979, z v = 1.021. The sample of haloes at z = 
consists of 36500 haloes each with mass M = W^h^MQ. 
The number of haloes in the range 1.5 - 2.5 x 10 13 /i _1 M Q 
at z d is N p = 19241 and N d = 15465. In the range 
4.5 - 5.5 x 10 13 /i _1 M Q there are N p = 7587 and N d = 4961 
haloes. Finally, in the range 6.5 — 7.5 x 10 13 /i _1 Mq there 
are N p = 5803 and N d = 3330 haloes. The resolution mass, 
Mmi n used is 5 x 10 11 hT 1 Mq. For z = 1 long dashes show 
the formula of SBBW08, dots the formula of FM08 and small 
dashes the prediction of the binary spherical model given by 
Eq. 6. In the following we examine the role of some -of the 



Figure 2. Thin solid line, dashed-dot-dot line and thick solid 
line show the mean merger rate at z f» 1 for haloes with masses 
in the range M 1 = 1.5 - 2.5 X 10 13 h' 1 M Q , M 2 = 4.5 - 5.5 X 
lO 13 fc -1 M0 and M 3 = 6.5-7.5 X 1O 13 /i~ 1 M , respectively. The 
merger tree used started at z = with D = 0.05 for a sample 
of 36500 haloes each of mass lO 14 /i _1 Af . After 22 time steps 
the redshifts are z d = 0.979, z p = 1.021. The resolution mass 
is M m i„ = 5 X lO 11 ^ - 1 Mq. Long dashes, big dots and small 
dashes are the results from the formulae of SBBW08, FM08 and 
the binary spherical model for M d = 5 X 10 13 /i — 1 Mq at z = 1. 



large number- of the parameters involved in the construction 
of merger trees. These are: 

(i) the resolution mass M m in 

(ii) the step in redshift, defined by the parameter D 

(iii) the value of as and 

(iv) the number of realizations, that is the number of 
haloes at 2 = 0. 

Resolution mass M m in is a crucial parameter. A sample 
of present day haloes of mass Mo is analyzed to smaller 
and smaller haloes as the redshift becomes higher. Let M d 
be the mass of a descendant halo at some high redshift, 
Mi its largest progenitor and M2 another progenitor. 
Obviously, Mi < M d - M m in and M 2 > M m4n . Thus 
£ = iff ^ M d M Mm»n = 5™" and consequently £ ml „ and 
Mmin are related by M min = M d (l + C m L) _1 - Since at high 
redshifts M d is significantly smaller than Mq, the condition 
for the merger rate curves to extend to the left up to values 
as small Mm^n < M ( 1 + CD" 1 - We found 

that for decreasing M m i„ the value £, at which the merger 
rates curves show their rapid fall, moves to the left. Thus, 
this rapid fall seen in Figs 1 and 2 is clearly a matter of 
resolution. 

On the other hand, the results do not seem to be sensitive 
to the values of the step in redshift z. This step depends 
on the parameter D described in section 2.1. We examined 
cases, for D=0.025, D=0.05 and D=0.1. We found that the 
results are the same. 

Fig. 3 shows the role of the above two parameters. It presents 
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Figure 3. Merger rates for descendant haloes with masses in 
the range 1.5 — 2.5 X 10 13 /i _1 Mq, for an initial sample of 10000 
haloes each with mass I0 14 h~ 1 Mq after 72 steps (thin solid line), 
after 144 steps (thin solid line with large dots) and after 36 steps 
(line with thin and small dashes). These three lines correspond to 
about the same redshift (z=3), and to the same resolution mass 
M m in = 5 X 10 11 h^ 1 Mq. Thick solid line is the prediction for 
smaller minimum mass, M mi „ = 5 x 10 w h- 1 M Q . Big dots arc 
the predictions of FM08, while large thick dashes, the predictions 
of SBBW08. 

merger rates, for descendant haloes with masses in the 
range 1.5-2.5 x 10 13 h~ 1 M Q , for an initial sample of 10000 
haloes, each with mass 10 4 /i _1 Mq, after 72 steps (thin 
solid line), after 144 steps (thin solid line with large dots) 
and after 36 steps (line with thin and small dashes). All 
lines correspond to about the same redshift (z=3), since we 
used different values for the time step parameter, D — 0.05, 
D = 0.025 and D — 0.1, respectively. The corresponding 
redshifts are (z d = 2.971, z p = 3.011), (z d = 3.005,%, = 
3.025), (z d = 2.921, z p = 3.000), respectively. We see that 
the resulting merger rates are similar, although the values 
of dz used in Eq. 3 differ, dz w 0.02, 0.04, 0.08. 
In the same figure the thick solid line is the prediction for 
a resolution mass an order of magnitude smaller, namely 
M m i n — 5 x 10 10 /i~ 1 Mq and shows the above described 
role of M m in- Finally, big dots and thick dashes are the 
predictions of FM08 and SBBW08, respectively, for values 
of the parameters given above. All calculations in Fig. 3 are 
for 10000 present day haloes and for as = 1. The above 
described role of the two parameters is the same for the 
range of present day masses, W 10 Mq < M a < 10 15 Mq, 
we examined. We note that small values of M m in lead 
to large numbers of haloes at the past. For example, 
10000 present day haloes with mass 1O 15 M have for 
Mmin = 5 x 10 12 Mq, about 1.2 x 10 5 progenitors at z w 3. 
For the minimum value of M m in that we have used, that 
is M m i n = 5 x 10 10 Mq and corresponds to a resolution 
1 : 20000, at the same redshift, the number of progenitors 
is larger than 1.5 x 10 6 . Thus, tree construction becomes a 
computing time consuming procedure. 



Figure 4. Merger rates for descendant haloes with masses in the 
range 1.5 — 2.5 x 10 13 h~ 1 MQ at z = 1. Dots are the predic- 
tions for the power spectrum normalized to erg = 1, while dashes 
correspond to as = 0.9. Solid line shows a law oc £ -1 - 52 for com- 
parison. Deltas correspond to the model of FM08 while solid line 
with black squares to the model of SBBW08. 



Major mergers seem to play an important role in the 
formation of dark matter haloes. N-body simulations show 
that haloes which experienced a recent major merger 
event, appear to have lower concentrations and steeper 
inner density profiles (e.g. Ascasibar et al. 2003, Tasitsiomi 
et al. 2004) as well as larger values of spin parameters 
(Gardner 2001, Peirani et al. 2004). This last result is also 
supported by semi-numerical results (e.g. Vitvitska et al. 
2002, Hiotelis 2008). Thus, it is interesting to study in more 
detail the behavior of merger rates curves at large values 
of £. This demands smooth curves and consequently large 
number of haloes. So, in Fig. 4 we present an example of 
the role of the last two of the parameters we examined. 
These parameters are the number of present day haloes 
and the value of as. All curves of this figure correspond to 
2 = 1. Dots are the predictions for a set of 36500 present 
day haloes each with mass Mo = 10 14 h~ 1 Mq for a power 
spectrum with as = 1 and for descendant haloes in the 
range 1.5 x IO^/i^Mq - 2.5 x 10 13 /i _1 M Q . At z = 1, 
there are N p — 19241 and N d = 15456 haloes. Dashes are 
the predictions for a set of 1.5 x 10 present day haloes 
of the same mass for as = 0.9 and for the above range 
of mass of the descendant haloes. At 2 = 1, there are 
N p = 104873 and N d = 70276 haloes. Solid line shows a law 
oc £ -152 for comparison. Deltas are the predictions of the 
fitting formula of FM08 and finally, the solid line with big 
dots shows the predictions of SBBW08. Fig. 5 is similar to 
Fig.4 but for 2 = 3 and for descendant haloes in the range 
5 x 10 12 h~ 1 M Q - W 13 h- 1 M Q . There are N p = 26610 and 
N d = 24540 for the tree corresponding to the dashed line 
while there are N p = 129074 and N d = 109678 for the 
tree corresponding to dots. It is clear from the above two 
figures, that - at the level of accuracy of the calculations 
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Figure 5. As in Fig. 4 but for z=3. 

of this paper - no differences can be detected due to the 
difference between 0.9 and 1.0 of the values of erg. Larger 
number of present day haloes leads, obviously, to smoother 
curves but they cause no difference in the overall shape. 



The construction of reliable merger-trees is a subject under 
current investigation. These algorithms usually have funda- 
mental problems regarding mass conservation, accurate rep- 
resentation of the distribution of progenitors etc. As regards 
N-body simulations, in addition to their resolution prob- 
lems, it is characteristic that some of their results depend 
on the techniques used for their derivation. For example, 
Bett et al. (2007) showed that the values of spin parame- 
ters of haloes and their behavior as a function of mass de- 
pends crucially on the halo-finding algorithm. However, the 
study of physical parameters -as for example a more accu- 
rate power spectrum- seems not to be proper at this stage. 
Probably, the results are more sensitive to the method used 
for their analysis than to physical parameters and thus, any 
detailed description between the results of N-body and those 
of merger-trees may not be very useful. However, there are 
characteristic trends that show interesting agreement. 
Further improvements of tree construction algorithms as 
well of the quality of N-body simulations could help to un- 
derstand better the physical picture during the process of 
the formation of dark matter haloes. 
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4 DISCUSSION 

The above study of merger rates that are predicted by 
merger trees based on the extended PS theory using the 
ellipsoidal collapse model, leads to the following conclusions 
that hold for haloes in the range of mass 10 12 /i _1 Mq — 
10 15 /i -1 Mq and for redshift < z < 3. 

(i) Merger rates depend on the mass of the descendant 
halo as ~ M™ with n= 0.2. This dependence is in practice 
indistinguishable from the value n = 0.15 proposed from 
SBBW08 but is far from the value of n = 0.083 proposed by 
FM08. 

(ii) We found that merger rates depend on the redshift 
through the quantity [dS c (z)/dz] 1 , where I — 1.1. This value 
of the exponent I is close to the value predicted by the binary 
spherical case (see Eq.6) and between the values 0.371 and 
2 proposed by FM08 and SBBW08, respectively. 

(iii) The results of merger trees are in better agreement 
(for the above values of n and I) with the predictions of for- 
mulae of FM08 and SBBW08 than with the predictions of 
binary spherical model given by Eq.6. The binary spheri- 
cal model underestimates merger rates for small values of £ 
while merger trees, for proper resolution, give results much 
closer to those of N-body simulations. 

(iv) For large values of £, merger rate curves scale approx- 
imately as w £ . 

(v) We examined steps in redshift from ~ 0.02 to w 0.08 
and we found no differences in the results. 

(vi) Smaller values of the resolution mass M m i n give a 
better agreement with the above fitting formulae for small 
values of £. 
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